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ABSTRACT 

We estimate the spectral properties of the stochastic background of gravitational ra- 
diation emitted by a cosmological population of hot, young, rapidly rotating neutron 
stars. Their formation rate as a function of redshift is deduced from an observation- 
based determination of the star formation history in the Universe, and the gravita- 
tional energy is assumed to be radiated during the spin-down phase associated to the 
newly discovered r-mode instability. We calculate the overall signal produced by the 
ensemble of such neutron stars, assuming various cosmological backgrounds. We find 
that the spectral strain amplitude has a maximum m (2 — 4) x 10~ 26 Hz~ 1//2 , at 
frequencies « (30 — 60) Hz, while the corresponding closure density, h?Q,GWi has a 
maximum amplitude plateau of « (2.2 — 3.3) x 10~ 8 in the frequency range (500 — 1700) 
Hz. We compare our results with a preliminary analysis done by Owen et al. (1998), 
and discuss the detectability of this background. 
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1 INTRODUCTION 

In the last two years, a series of investigations on the per- 
turbations of rotating relativistic stars have revealed the ex- 
istence of a class of r-modes that are unstable due to the 
emission of gravitational radiation (Andersson 1998; Fried- 
man & Morsink 1998; Lindblom, Owen & Morsink 1998). 
This gravitationally driven instability is of considerable im- 
portance because it has a number of interesting astrophysi- 
cal implications. For instance, it determines a spin down of 
newly born neutron stars to periods closer to the initial peri- 
ods inferred for ms pulsars. Moreover, as the star spins down 
an energy equivalent to w 1% of a solar mass is emitted in 
gravitational waves making the process a very interesting 
source for gravitational-wave detection (Owen et al. 1998; 
Andersson, Kokkotas & Schutz 1998). 

In this paper we study the stochastic background of 
gravitational waves contributed by a cosmological popula- 
tion of hot, young, rapidly rotating neutron stars with a 
small initial excitation in the I — 2 r-mode. A rough esti- 
mate of the spectrum of this background radiation has re- 
cently been given by Owen et al. (1998), assuming a comov- 
ing number density of neutron star births which is constant 
in the range < z < 4 and zero at earlier times. 

Here we present a more complete analysis, which derives 



the neutron star birth rate as a function of redshift from an 
observation-based determination of the star formation den- 
sity evolution (Madau, Pozzetti & Dickinson 1997), within 
three different cosmological background models. Thus, the 
present investigation relies on a source rate which correctly 
accounts for the observed strong luminosity evolution of star 
forming galaxies between z = and z — 1 — 2 (Lilly et al. 
1998; Connolly et al. 1997). With a similar approach, the 
gravitational background emitted by an ensemble of stellar 
mass black holes has been recently computed by Ferrari, 
Matarrese & Schneider (1998). 

The plan of the paper is as follows. In Section 2 
we calculate the universal rate of neutron star formation; 
in Section 3 we briefly introduce the adopted model of 
gravitational- wave emission by an unstable rapidly rotating 
neutron star, closely following the recent analysis by Owen 
et al. (1998); in Section 4 we estimate the overall stochas- 
tic background produced by the ensemble of such sources 
throughout the Universe and calculate the total duty cy- 
cle of the process; Section 5 contains a brief analysis of the 
detectability of such a background, while conclusions are 
drawn in Section 6. 
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2 THE RATE OF NEUTRON STAR 
FORMATION IN THE UNIVERSE 

Type II and Type Ibc supernovae are fuelled by the gravi- 
tational contraction of their core (e.g. Ruiz-Lapuente 1997). 
Numerical studies have shown that single stars with masses 
between 8 and 20Mq form iron cores with masses near 
the Chandrasekhar limit (e.g. Timmes, Woosley & Weaver 
1996). Thus, there is universal agreement that these stars 
leave neutron star remnants. The fate of stars with masses 
greater than 20M© is less certain and depends on the iron 
core masses formed and on the effect of fallback during the 
explosion (e.g. Woosley & Weaver 1995; Woosley & Timmes 
1996). 

For the calculations to follow, we assume that stars with 
masses between 8 and 25Mq give origin to neutron stars. 
However, we will also investigate the effect of an upper mass 
cutoff of ~ 20M Q . 

The number of neutron stars formed per unit time 
within the comoving volume out to redshift z is given by 
(Ferrari, Matarrese & Schneider 1998) 

r dv r 25M ® 

R NS (z) = / p.(z') -^7 dz' / $(M) dM. (1) 

JO dZ J8M e 

Here p*(z) is the star formation rate density (mass of gas 
that goes into stars per unit time and unit comoving vol- 
ume), <I>(M) is the initial mass function (IMF), and dV/dz 
is the comoving volume element. 

Following Madau, Pozzetti & Dickinson (1997), we 
model the star formation evolution from the recent collec- 
tion of UV-optical observations of star forming galaxies out 
to redshifts ~ 4 — 5, by assuming a Salpeter IMF, 

$(M) oc M- {1+x) , with x = 1.35 , (2) 

normalized through the relation 

/■125M 

/ M $(M) dM = 1. (3) 
Jo.iM 

The star formation history is plotted in Figure 1 for a flat 
cosmology with vanishing cosmological constant and h = 0.5 
(here h is the Hubble constant, Ho, in units of 100 km s _1 
Mpc -1 ). The star formation rate rises sharply from its lo- 
cal value to a peak at z ~ 1.5 to fall again out to z ~ 4. 
This behaviour indicates that the bulk of the stellar popu- 
lation was assembled in the redshift range 1 ~ z ~ 2. This is 
consistent with the indications of QSO absorption lines and 
metallic clouds observations (Lanzetta, Yahil & Fernandez- 
Soto 1996; Pei & Fall 1995) and with the predictions of semi- 
analytic models for a broad class of hierarchical clustering 
cosmologies (Baugh et al. 1998). A thorough review of the 
current state of observations of the star formation history 
of field galaxies is given in Ellis (1997). A summary of the 
relevant aspects for the present analysis is given in Ferrari, 
Matarrese & Schneider (1998). 

The predictions for the rate of core-collapse supernovae 
are in good agreement with the observed local values (Cap- 
pellaro et al. 1997). At z > 1, the detection of Type II events 
must await future experiments, such as the Next Generation 
Space Telescope (Madau, Delia Valle & Panagia 1998). 

Here we shall only stress that one major source of un- 
certainty in the interpretation of the data points in Figure 
1 is the effect of dust extinction. Even a small amount of 




Redshift 

Figure 1. Evolution of the SFR density (M Q yr- 1 Mpc -3 ) for 
Q = 1, (A = 0), h = 0.5 and a Salpeter IMF. The data points 
are taken from Gallego et al. (1995) (empty dot), Treyer et al. 
(1998) (filled triangle), Lilly et al. (1996) (filled dots), Connolly 
et al. (1997) (empty squares) and the HDF (Madau et al. 1996) 
(filled squares). The solid line represents the SFR density evo- 
lution assumed in the present analysis. The fit has been kindly 
provided by P. Madau. 



dust in primeval galaxies can significantly attenuate the UV 
luminosity and re-radiate it in the far IR thus leading to a 
potentially serious underestimate of the star formation ac- 
tivity. Dust correction factors are still very uncertain mainly 
because of the unknown shape of the dust UV-extinction law 
(Pettini et al. 1997). The model for the global SFR evolu- 
tion presented in Figure 1 assumes an extinction law sim- 
ilar to that which applies to stars in the Small Magellanic 
Cloud (SMC) and an amount of dust which results in an up- 
wards correction of the comoving UV luminosity by a factor 
1.4 at 2800A and 2.1 at 1500A (Madau, Pozzetti & Dickin- 
son 1997). Recent near infrared observations of a galaxy at 
z — 4.92 has given evidence that absorption by dust is im- 
portant up to redshifts of ~ 5 (Soifer et al. 1998) . Similarly, 
analysis of the Infrared Space Observatory (ISO) data have 
shown significant departure from the UV-optically derived 
SFR evolution at z > 2 (Rowan- Robinson et al. 1997). 

This may seem a serious caveat for the present investi- 
gation. However, this is not necessarily the case because the 
extra-galactic astrophysical gravitational- wave backgrounds 
prove to be mainly contributed by low-to-intermediate red- 
shift sources. High redshift sources play a negligible role 
because the energy flux emitted by each source decreases 
as the inverse of the squared luminosity distance. On the 
other hand, if an upwards correction of the UV-optically 
derived SFR density is required at z w 1.5, as suggested 
by a recent detection of the infrared background by the Dif- 
fuse Infrared Background Experiment (DIRBE) and the Far 
Infrared Absolute Spectrometer (FIRAS) on board of the 
Cosmic Background Explorer (COBE) satellite (Dwek et al. 
1998), this would result in an upwards correction of roughly 
the same amplitude on our gravitational-wave background 
predictions. 

We have evaluated equation (|l|) for three sensible cos- 
mological models (e.g. Baugh et al. 1998), namely a flat ge- 
ometry with vanishing cosmological constant, an open back- 
ground and a flat low-density model (see Table 1). For a gen- 
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Table 1. Cosmological parameters adopted in the three back- 
grounds considered. 

Model Q M n A h 



A 


1 





0.5 


B 


0.3 


0.7 


0.(5 


C 


0.4 





0.6 



30 




1 2 3 4 5 

Redshift 



Table 2. Values of the total rate of neutron stars in different 
cosmological scenarios and for two values of the upper progenitor 
mass cutoff (20Af Q , 25M Q ). 

Model M up /Mq Rns I (eventspers) 



A 25 19.5 

A 20 17.6 

B 25 28.7 

B 20 25.9 

C 25 21.7 

C 20 19.6 



sive ones which, in turn, are those that at the end of their 
evolution give origin to core-collapse supernovae. 

The total number of supernova explosions per unit time 
leaving behind a neutron star is given in Table 1 for the con- 
sidered cosmological models, and for two values of the upper 
mass cutoff, M = 25Mq and M = 20Mq. As a consequece 
of the largest number of stars having masses < 25Mq, the 
total rate of core-collapses to neutron stars is a factor ~ 4 
larger than that of core-collapses to black holes considered 
in Ferrari, Matarrese & Schneider (1998). 



Figure 2. Evolution of the rate of neutron star formation for 
three different cosmological scenarios. 

eral background, the comoving volume element is related to 
z through, 

dV = 47t(^J r 2 T{n M ,n A ,z)dz, (4) 
where Ho is the Hubble constant and 

F(n M , n A , z) = 1 . (5) 

y/(l + Zf (1 + Sl M Z) - Z (2 + Z) Q A 

In eq. (Q) r is the comoving distance, 

r=— ^=s(VtT K f [(l + z'f(l + n M z')- (6) 
nov"K Jo 

z {2 + z')Q A ]~ 1/2 dz') 

and the function S has the following expression (Kim et al. 
1997) 

J]m + ^a>1 S(x) =sin(z), £l K = 1 - CIm + A , (7) 
Qm + f^A < 1 S(x) = sinh(a;) fi K = Qm + Q A — 1, 

Om+^a = 1 S(x)=x n K = i. 

Consistently, the SFR history shown in Figure 1 needs to be 
properly rescaled in order to account for the cosmological 
background dependence of both the luminosity distance and 
the comoving volume within each redshift bin. 

By using equation (^) we have evaluated the function 
Rns(z) defined in eq. ([j]) for the three cosmological models, 
and the results are shown in Figure 2. 

It should be mentioned that the evaluation of the rate 
of neutron star formation, based on the UV-optically esti- 
mated star formation rate evolution p*(z), and on an as- 
sumed universal IMF, is largely independent of the specific 
IMF adopted (Madau 1998). In fact, the stars which mainly 
contribute to the UV-optical luminosity are the most mas- 



3 GRAVITATIONAL EMISSION FROM 

RAPIDLY ROTATING NEUTRON STARS 

In this section we shall briefly describe how the timescale of 
the spindown phase and the energy spectrum of the gravita- 
tional radiation emitted during the r-mode instability have 
been estimated (Lindblom, Owen & Morsink 1998; Owen et 
al. 1998; Andersson, Kokkotas & Schutz 1998; Andersson, 
Kokkotas, Stergioulas 1998). These quantities will be used 
to compute, respectively, the duty cycle and the spectral 
energy density of the stochastic background. 

The excitation of an r-mode starts as a small pertur- 
bation of the velocity field. As the amplitude of the mode 
grows, non-linear hydrodynamic effects become important 
and eventually dominate the dynamics. Under the assump- 
tion that the star is uniformly rotating with angular velocity 
Q, and by parametrizing the r-mode by its amplitude a, the 
star can be considered as a system having only two degrees 
of freedom, Q and a. In this simple model, the total angular 
momentum for the i — 2 r-mode can be shown to be 

J(Q,a) = (J- | Ja) Q.MR 2 , (8) 

where 

i =zm?j Q pr4dr and J = WW 2 j pr " dr ' (9) 

and the rotational energy is 

E^J{Q,a)Q.. (10) 

Following Owen et al. (1998) we shall consider a New- 
tonian model of neutron star composed of a fluid with a 
polytropic equation of state, p = k p 2 , with k chosen so that 
M = 1.4M Q and R = 12.53 km. For this model, 7 = 0.261 
and J = 1.635 x 10~ 2 . It should be noted that I and J 
depend only on the chosen polytropic index. 

The evolution of the angular momentum of the star is 
determined by the emission of gravitational waves, which 



© 1998 RAS, MNRAS 000, §-| 



4 V. Ferrari, S. Matarrese and R. Schneider 



couple to the r-modes through the current multipoles, pri- 
marily that with I = m = 2. For this mode, the frequency 
of the emitted gravitational radiation is v — (2/3n) fl. The 
evolution of fi and a during the phase in which a is small, 
can be determined from the standard multipole expression 
for angular momentum loss, and from the energy loss due 
to the gravitational emission and to the dissipative effects 
induced by the bulk and shear viscosity. In this phase f2 
is nearly constant and a grows exponentially. After a short 
time, that Owen et al. estimate to be of order 500 s, the am- 
plitude of the mode becomes close to unity and non-linear 
effects saturate and halt a further growth of the mode. This 
phase lasts for approximately 1 yr, during which the star 
loses angular momentum radiating approximately 2/3 of its 
initial rotational energy in gravitational waves, up to a point 
where the angular velocity and the temperature are so low 
that the r-mode ceases to be unstable. Viscous forces and 
gravitational radiation damp out the energy left in the mode, 
and the star slowly reaches its final equilibrium configura- 
tion. 

Since gravitational radiation is emitted over such a long 
time interval, the background contributed by neutron stars 
is continuous. In fact, the total duty cycle for this process is 



Vmin ~ (2/3tt)Q, c 



(16) 



10 J 



(11) 



D= / dR NS (z)T NS (l + z) 
Jo 

It has been suggested that if the proton fraction in the star is 
sufficiently large to make the direct URCA process possible, 
the star would cool in a much shorter time, rjvs ~ 20 sec. 
However, if this were the case, the duty cycle would still be 
3> 1 (D ~ 10 2 ) and the resulting background would still be 
continuous. 

It is known that the angular velocity of a rotating star 
cannot exceed the Kepler frequency fijc at which mass shed- 
ding at the stellar equator makes the star unstable. A rea- 
sonable approximation for this limiting frequency is 



C 



(12) 



M / 10 km\ 3 
M© V R ) ' 

where C w 7.8 x 10 3 sec -1 (Friedmann, Ipser & Parker 1989), 
and M and R refer to the mass and radius of the corre- 
sponding non-rotating star. Following Owen et al., we shall 
assume that neutron stars are born with angular velocity 
close to the maximum value Qk, and have a rotational ki- 
netic energy 



E K « J {fix; a) U K 



(I ~^Ja 2 ) M R 2 Sf K . 



(13) 



In the non-linear saturation phase, where most of the grav- 
itational radiation is emitted, awl. Since in this simple 
model approximately 2/3 of the stars rotational energy is 
radiated in gravitational waves, the expression of the energy 
spectrum can be approximated as follows 

dEaw 4 



dv 
where 



(14) 



v-max = (2/37r)fijr- (15) 

The r-mode instability ceases to be effective for some value 
of the rotational frequency Q c , which corresponds to the 
final spin period, and which is associated to a gravitational 
emission frequency 



Q c can be determined by solving the equation l/r(fi c ) = 0, 
where r is the total dissipation timescale which can be de- 
composed as a sum of the damping times associated to the 
gravitational emission, to the shear and to the bulk viscos- 
ity. t(Q c ) is clearly a function of the temperature of the star, 
and it has been shown that the r-mode instability operates 
only in hot neutron stars (T > 10 9 K). Below 10 9 K, super- 
fluidity and other non-perfect fluid effects become important 
and the damping due to viscosity dominates with respect to 
the destabilizing effect of the gravitational radiation. To be 
consistent with the analysis by Owen et al. (1998), we shall 
adopt a value of f2 c m 566 Hz for the final angular velocity 
which corresponds to a final spin period of w 11 ms and to 
Vmin ~ 120 Hz. 

The qualitative picture that arises from this simple 
model is believed to be sufficiently reliable, even though 
various uncertainties and approximations might affect the 
quantitative results for the initial rotation of the star after 
collapse, for the spindown timescales as well as for the final 
rotation period ( Andersson, Kokkotas & Schutz 1998) . Thus, 
it will be interesting to evaluate the effects that different val- 
ues of v m ax, Ek and v m in have on our main conclusions. 



4 SPECTRAL PROPERTIES OF THE 
STOCHASTIC BACKGROUND 

In order to evaluate the spectral energy density of the 
stochastic background produced by the spindown radiation 
from newly born neutron stars, we need to know how the 
flux emitted by a single source at redshift z would be ob- 
served today, and convolve it with the differential rate of 
the sources dRNs(z) as in Ferrari, Matarrese & Schneider 
(1998). 

The average energy flux per unit frequency emitted by 
a single neutron star at cosmological redshift z is 



1 f dEgw Y 
4TYd L (z) 2 \ dv ) 



(17) 



where v obs is the frequency at which the radiation would 
presently be observed. The energy flux f(v° s ) can be evalu- 



ated by making use of eq. ( |14[ ) suitably redshifted as follows 

obs 

(18) 



£/ otss Ek 



1 



3tt (v^S x y (l + z)d L (z) 2 



where 
v em 
l + z 



em 

obs ^max 

^v «S — — , 

1 + 2 



(19) 



and di,(z) = (1 + z) r is the luminosity distance, which de- 
pends on the cosmological model as indicated in eq. (Q). 
Thus, the spectral energy density of the produced back- 
ground is 



dEaw 
dtdSdv 



f(v) dRNs(z 



where the differential rate of sources is given by 

/•25M Q 



dRNs(z) = p»(z) — 



$(M) dM 



8A/ :: . 



dz. 



(20) 



(21) 
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Figure 3. The spectral energy density (dE / dSdtdv) is plotted 
as a function of the observational frequency for the cosmological 
background corresponding to model A (see Table 1). 
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Figure 4. The spectral strain amplitude S h corresponding to 
the spectral energy density given in Figure 3 is plotted as a func- 
tion of the observational frequency. 



[see eq. (gj)]. From the spectral energy density the corre- 
sponding values of the closure energy density of gravitational 
waves 

v bs dE 



c 3 p cr dtdSdu' 



(22) 



where p cr = 3Hq / 8tv G, and of the spectral strain ampli- 
tude 



r— (2G 1 \ 1/2 f dE y/ 2 



(23) 



can be derived. 

In Figures 3 and 4 we plot the spectral energy density 
and the strain amplitude for the flat background with zero 
cosmological constant corresponding to model A (see Table 
1). The effect of a varying cosmological background is neg- 
ligible on both the spectrum and the strain amplitude. In 
fact, the amplification of the rate at high redshifts shown 
in Figure 2 for the models B and C, is mostly suppressed 
by the inverse squared luminosity distance dependence of 
the single source spectrum [see eq. (EsJ)] for the same mod- 
els. Conversely, the closure density exhibits a more evident 
dependence on the Hubble constant, as shown in Figure 5, 
where we have considered only model A and B because the 
dependence on the other cosmological parameters (namely 
Qm and Qa) is negligible. 

The estimate of the intensity of the stochastic back- 
ground also depends on the chosen values of f m in and v m ax- 
We have repeated our calculations for two values of v m in, 
chosen so as to cover the range of pulsar rotational periods, 
Pi n , inferred from observations 



-P-- 1 



(24) 



with Pin between 9 and 22.8 ms, where the upper limit is 
the maximum rotational period at which the r-modes are 
unstable (Andersson, Kokkotas & Schutz 1998). Since the 
strain amplitude is sensitive mainly to the low frequency re- 
gion of the spectrum, due to the factor v~ b 2 s in eq. (|23|), a 
decrease of v m in determines an increase of the peak am- 
plitude and shifts its location towards lower frequencies, as 
shown in Figure 6. 



o 
C! 




600 800 
Frequency (Hz) 

Figure 5. The function Ogw corresponding to the spectral en- 
ergy density plotted in Figure 3. Since the main effect on the grav- 
itational background is due to the change in the Hubble constant, 
we have plotted only two curves corresponding to an h = 0.5 fiat 
model without cosmological constant and to an h = 0.6 flat model 
with a cosmological constant (see Table 1). 



The value of i/ m „, depends on Qk which, in turn, de- 
pends on the mass and radius of the star, as shown in eqs. 
( [iB] ) and (|l2j). The plots shown in Figures 3, 4 and 5 were 
obtained by assuming M — 1.4Mq and R — 12.53 km. In 
reality, the possible radii for rotating neutron stars range 
from 10 to 15 km, therefore, since the dependence of Qk 
on the mass is less significant than that on the radius, we 
have computed the strain amplitude and the closure den- 
sity by considering M — IAMq, and changing the radius to 
R — 10 km and to R — 15 km. The effect of reducing v max , 
i.e. of increasing the radius, is that of narrowing the interval 
of frequencies over which each source contributes. Since the 
energy flux of a single event does not depend on £Ik, [see 
eq. (US)], this will result in an increase of the amplitude of 
the peak of the spectral energy density. This effect is am- 
plified in the strain amplitude, that is more sensitive to low 
frequencies. However, it should be stressed that the increase 
is smaller than a factor 2. Conversely, the closure density 
is more sensitive to the high frequency contributions, and 
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Figure 6. The spectral strain amplitude S, corresponding to a 
flat model is plotted as a function of the observational frequency 
for three different values of v m i n (see text). 
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Figure 7. The spectral strain amplitude S h corresponding to a 
flat model is plotted as a function of the observational frequency 
for three different values of the radius and, thus, of the initial 
rotational period and of the total rotational energy (see text). 



therefore the effect is the opposite. The results are shown in 
Figures 7 and 8. 



5 DETECTABILITY 

The value of the optimal signal-to-noise ratio (SNR) for re- 
shaped interferometers can be computed according to the 
formula (Allen 1996) 



SNR" 



V 50tt 4 I 



T 



dv 



7 2 M!iGw(^) 



v 6 S 



(i) 



(25) 



where are the spectral noise densities of two detec- 

tors operating in coincidence, and is the overlap func- 

tion which takes into account the difference in their location 
and orientation. We shall calculate SNR for several pairs 
of detectors, and for £Igw computed for a fiat cosmology, 
c mm = 120 Hz, and v max = 1396 Hz. We shall assume an 
integration time of T = 1 yr. 
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Figure 8. The function £Igw corresponding to the spectral strain 
amplitude plotted in Figure 7. 



For interferometric antennas, we have computed 
from the expected sensitivity curves given for VIRGO (Bec- 
caria et al. 1996) and for ADVANCED LIGO (Flanagan 
& Hughes 1997). The sensitivity curve for GEO600 has 



been extract ed from the Web site http://www.geo600.uni- 



hannover.de 



For VIRGO and GEO600 we find SNR = 3.7 x 10~ 3 , 
assuming optimal orientation, and SNR — 5.0 x 10 -3 if 
GEO600 is narrowbanded around 100 Hz. For two AD- 
VANCED LIGO detectors SNR = 0.84. Since the corre- 
lation of two interferometers at a distance d starts to be 
very poor around a frequency ~ (70 Hz)(3000 km/d) (Allen 
& Romano 1998), two nearby detectors respond better to 
higher frequencies, where our predicted signal has a larger 
amplitude (see Figure 8). Thus, the signal to noise ratio 
rises to 8.5 if we consider two nearby ADVANCED LIGO 
(or VIRGO) interferometers at a distance of ~ 300 km. 

We have similarly investigated how interferometric and 
resonant detectors operating in coincidence might respond 
to our predicted background. The correlation of VIRGO 
and NAUTILUS, situated 265 km apart, gives SNR = 
3.0 x 10 -4 . Assuming that a similar correlation is done with a 
truncated icosahedral antenna (TIGA) (Merkowitz & John- 
son 1995) and an ADVANCED LIGO interferometer located 
at the same site, the resulting SNR is ~ 1.3 x 10 -3 , whereas 
it is SNR « 1.7x 10" 3 if the interferometer is ADVANCED 
VIRGO. 

Thus, the detection of this background requires two 
nearby interferometers with 'advanced' sensitivities, as sim- 
ilarly noted by Owen et al. (1998). 



6 CONCLUSIONS 

In this paper we have evaluated the gravitational back- 
ground produced by a cosmological population of hot, young 
and rapidly rotating neutron stars, emitting gravitational 
radiation during the spin-down phase associated to the r- 
mode instability. The dependence of this background on the 
cosmological model and on the parameters that enter into 
the determination of the spectrum of each single event, have 
been discussed. Our predictions for the spectrum, spectral 
strain amplitude and closure density, appear to be rather 
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Figure 9. A comparison with the closure density, h 2 £Igw> P re ~ 
dieted by Owen et al. (1998) is plotted as a function of the obser- 
vational frequency for the cosmological background corresponding 
to model A. 



insensitive to the assumed cosmological parameters. This is 
because the increase of the neutron star birth rate we find 
in an open or a flat low-density cosmology is compensated 
by the increase of the luminosity distance as a function of 
redshift, which suppresses the contribution of the farthest 
sources. 

For the same reason, the uncertainties that still trouble 
the high redshift tail of the SFR density evolution due to 
dust obscuration, are not severe for our signal, as this is 
mainly contributed by low-to-intermediate redshift sources. 
Moreover, if dust obscuration produces an upward correction 
of the SFR density at intermediate redshifts (i.e. around the 
location of the maximum) then the amplitude of (Sh) 1 ^ 2 and 
of Qaw will be similarly amplified at ~ 30 Hz and ~ 600 
Hz respectively. 

We find that, within a reasonable range of values of the 
main parameters which characterize the energy spectrum of 
a single source (namely v m in, Vmax and EjS), the spectral 
strain amplitude has a maximum ranging within ~ (2 — 4) x 
10" 26 Hz~ 1/2 at frequencies « (30 - 60) Hz, whereas the 
closure density, h 2 Q.gw , is shown to have a long plateau 
extending from w 300 Hz up to « 1700 Hz with an amplitude 
of « (2.2 - 3.3) x 10 -8 . 

It is interesting to compare our results with those of 
a similar study done by Owen et al. (1998), who assume a 
constant comoving number density of neutron star births in 
the range < z < 4, and zero elsewhere. 

In Figure 9 we show the closure density predicted by the 
present analysis and that obtained by Owen et al. The two 
curves exhibit quite a different high frequency behaviour. 
Since the cosmological expansion shifts the frequencies to- 
wards smaller values, the high frequency tail corresponds to 
closer sources, thus the difference between the two curves 
can be traced back to the SFR density evolution between 
z = and z ~ 1 — 2. The peak at ~ 600 Hz exhibited 
by our curve, corresponds to z « 1.3, where we expect to 
have the maximum number of sources. Conversely, the max- 
imum amplitude predicted by Owen et al. is located around 
280 Hz, and corresponds to the energy emitted at maximum 
spin rate by a source at z ~ 4. Thus, their signal appears to 
be mostly contributed by distant sources near their maxi- 



mum rotational velocity, and this is due to the fact that they 
do not consider the luminosity distance damping, which is 
particularly effective for high redshift sources. It is impor- 
tant to stress that in the region where cross-correlation be- 
tween two LIGO interferometers can be accomplished, i.e. 
for v < 50 — 60 Hz, the signal is entirely emitted at z > 1, 
where it is important to embed the gravitational sources in 
the correct cosmological environment. 

Both in our analysis and in the study of Owen et al., all 
neutron stars have been assumed to be initially rotating with 
an angular velocity close to the Keplerian value, while there 
might be a fraction of stars born with lower initial spin rate. 
In particular, no spindown gravitational radiation would be 
emitted if the initial angular velocity of the star is too small 
for the instability to set in Spruit & Phinney (1998). How- 
ever, according to Andersson, Kokkotas & Schutz (1998) this 
is unlikely to be the generic situation, because the angular 
momentum of the degenerate core required to make the star 
rotating with velocity Qk is rather small. Moreover, there is 
clear evidence that some pulsars were born spinning rapidly: 
as an example, the observed spin period for the young X- 
ray pulsar in the supernova remnant N157B is 16 ms, and 
its inferred initial spin period is estimated to be « 6 — 9 ms 
(Marshall et al. 1998). 
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